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A Mars precision landing requires a landed footprint of no more than 100 meters. Obstacles to reducing the 
hnded footprint include 8 trajectory dispersions due to initial atmospheric entry conditions such as entry 
anele narachute deployment height, environment parameters such as wind, atmospheric density, parachute 
dep lo v ni en t ^dy nam lc ^unavoidable injection error or propagated error from launch, etc. Computational 
Intelligence (Cl) techniques such as Artificial Neural Nets and Particle Swarm Optimization have been 
shown to have great success with other control problems. The research period extended previous work on 
investigating applicability of the computational intelligent approaches. The focus of this investigation was 
on PaTc e IwTm Optimization and basic Neural Net architectures. The rese-r-h mvestrgating these issues 
was performed for the grant cycle from 5/15/01 to 5/15/02. Matlab 5.1 and 6.0 along with NASA s POST 
were the primary computational tools. 


Nomenclature 

Activation Function 
Artificial Neural Net 
Backpropagation 
Bank Angle 

Computational Intelligence 
Epoch 

Evolutionary Computation 

Firewire 

Generalization 

Genetic Algorithms 

Hyperspace 

Inverse Control 

Iris Data Set 

LM 

Mars Sample Return Mission 

MarsGram 

Matlab 

Multilayer 

Open Loop 

Particle Swarm Optimization 

POST 

Sigmoid 


Important Files Used 

Data & input: 

| bkb3guess.dat 
| chutestates.dat 
| fulliris.txt 
j prntblk.bak 
| prntblk.dat 
j pvstates.dat 
j bc4batch.mat 
j chutestatplot.mat 
| copy_of_bc 1871. mat 
| copy_of_lc.mat 

I copy_of_lc 1747 1 8. mat 

| copy_of_m200 1 ne wbase . mat 
| copy_of„rwplrps.mat 
j for ward, mat 
| for ward 3. mat 
j inverse, mat 
j IrisWtsLM.mat 
| lc4batch.mat 
j lcsaveas.mat 

| nettrain-76init-bcl7baseline.mat 

| nettrai n-7 6i nit-bc 1 7 basel ineMOD. mat 
| nettrain-nominit-bc 17baseline.mat 

| nettrain-nominit-bc 17baselineMOD.mat 
| nettrainer- weights. mat 
| nettrainer- wts-1 layer, mat 
| nettrainer2-wts- 1 layer-inertialBPX.mat 
| nettrainer2-wts- llayer-inertiallm.mat 
| nettrainer2- wts- 1 layer-inertialPSO.mat 
j pathdatl.mat 


Program & Control: 
chutestatplot.m 
demonormaiize.m 
f6.m 

learnga.m 

learnlm.m 

load_DATA_mat.m 

m2001viewsite.m 

MatPathData.m 

nettrainer2.m 

normalize. m 

p2check.m 

p3dbatchSPACE2DOO.m 

p3dbatchspace2000PSO.m 

pso.m 

testf6.m 

testnettrainer.m 

testpso.m 

tga2.m 

tlm2.m 

tpso2.m 

trainbpx.m 

trainga.m 

trainlm.m 

trainpso.m 

unwrapmat.m 

wrapmat.m 


Introduction 

The Mars Rover Sample Return Mission (MRSR) requirements make the development of an accurate 
descent control system challenging. A Mars precision landing requires a landed footprint of no more than 1 
kilometer. In contrast, the Viking mission had a landing footprint of over 100 kilometers (excellent for that 
particular mission requirement). The current state of the art’s footprint is still in the magnitude of 
kilometers. 

Obstacles to reducing the landed footprint include trajectory dispersions due to initial atmospheric entry 
conditions (entry angle, parachute deployment height, etc.), environment (wind, atmospheric density, etc.), 
parachute deployment dynamics, unavoidable injection error (propagated error from launch on), etc. 

In recent years, a number of computational techniques that lend themselves particularly well to control 
problems have become available. In broad categories they are Fuzzy Systems (FS), Artificial Neural 
Networks (ANN), and Evolutionary Computation (EC). Fuzzy Systems are based on fuzzy logic which is in 
turn based on the way our mind deals with incomplete and/or inaccurate information. Neural Nets are 
modeled after the spatial structure of the brain and allow ‘connectionist’ learning properties. Evolutionary 
computation paradigms such as Genetic Algorithms (GA) and Particle Swarm Optimization (PSO) are 
loosely modeled after biological evolution and optimization. These techniques have been shown 
individually to work very well in solving a large number of problems in linear and nonlinear system 
identification, modeling, and control. 

Until recently FS, ANN, and EC were totally separate fields with very little interaction. A growing number 
of researchers and practical engineers are discovering that a combination of two or more of these methods 
offers advantages that a single one lacks. 

Using an ANN in a control system we add fault tolerance, distributed (connectionist) representation 
properties, and the ability to learn optimal responses to new input. If we add an FS ‘shell’ we include high 
level rule/decision abilities as well as comparative reasoning. Combining techniques this way is called 
Computational Intelligence (Cl). 

This paper describes ongoing research into using Cl strategies for the MRSR mission as well as other Mars 
Lander missions. 

Problem Setup 

The problem examined is that of the re-entry of the proposed Mars 2001 lander. The lander specifications 
and landing site are all taken from the Mars 2001 conditions as related by Powell, Striepe, and Queen. The 
initial conditions are defined as within the cloud set developed by Powell. These initial conditions reflect a 
point of parachute deployment of the lander after it has completed the first leg of re-entry and consist of 
approximately 2000 possible starts within 10km dispersion. The 2000 points were reduced to 76 points by 
taking the convex hull of both the position and velocity subspaces. This results in a data set of initial 
conditions reflecting the extremes in handoff position and velocity. In addition, one initial point was 
computed from the averages of the complex hull and taken to be the nominal initial condition. These initial 
conditions represent the currently expected diversion of any Mars re-entry. 


Lander Scenarios 

The input deck setup to run with these initial conditions was a ballistic parachute deceleration followed by 
‘fire-wired’ steerable parachute descent. This scenario has been determined from previous work to have the 
best tradeoffs between practicality, stability, and performance. 

POST was used to target and optimize to desired end states. Because this input deck has been examined in 
the previous work only a brief overview is presented here. 

The lander simulation starts off with one of the 77 representative initial conditions. This includes all 
position and velocity data along with a time of day and basic weather model slightly modified from the 
MarsGRAM model. The modification is only in the wind strengths, using the Braun multipliers. 


The Post deck deploys the ballistic chute and at some controllable time it turns the ballistic parachute into a 
steerable parachute. Then several events follow where the controller is given the opportunity to change 
bank angle of the lander, all with a goal of minimizing distance to target. At a height 1000m above surface, 
where surface is taken to be at 2500m, the parachute is jettisoned and the lander performs a reverse gravity 
turn descent on thrusters to achieve desired end velocity. 

Controller Training 

This section describes motivation and previous research. 

For this deck the POST optimizer was turned off. The initial guesses for the controller (bank angles and 
times of occurrence) were generated in an offline data file to be fed to the deck on simulation start. For the 
nominal initial condition, which was the average of the extreme initial conditions, a set of controls was 
found by trial and error that would achieve fairly good results. These ‘nominal initial’ controls were then 
perturbed over a fairly large range for each initial condition, taking into account limits on bank angle, etc. 

For a situation where there was only one bank angle command at a controller specified time, we had 2 

controls: . 

1) time to turn parachute into steerable, also same as first bank angle command time 

2) bank angle command 

From previous work it was seen that a series of commanded bank angle turns yielded better results, so for a 
deck with 3 events we have 6 controls. Each set of these controls is perturbed parametrically and re-run on 
the same deck, i.e. for a deck with 2 control parameters the initial controls are run, then control 1 is 
changed slightly, deck is run again, then control 1 is changed again, and deck is run again. This is done 
until we get the whole range of possible (or at least reasonably expected) values for the control 1. Then 
control 2 is slightly changed and the whole parametric study of control 1 is re-done for the new control 2. 
The whole thing is repeated until we get a database reflecting the entire ranges of control 1 and 2. 

At the end of each run, all starting and ending trajectory data is saved. So for a deck with 6 controls and 
reasonably non-limiting control ranges, we have an immense amount of data to sift through. This data was 
used to train an artificial neural net as an inverse controller. The ending conditions along with initial 
conditions (inertial or relative position and velocity - both were tried) were the net’s inputs and the set of 
controls needed were the output. So for our 6 control deck we have a net with 12 inputs and 6 outputs. 

Not all of the collected data was used for training the net, about 20% was held back for testing after training 
was done. The trick to training is to make sure the data is well prepared. 

To this end a variety of things were tried. First, all the training data was normalized. This is a standard 
practice when using the traditional training algorithms. This required developing an algorithm that would 
normalize whatever data was given to it and save the normalizing factors so future unknown data could also 
be normalized in exactly the same way. This was coded as a matlab function where you could choose 
several types of normalization of scalars, vectors, or matrices. Also the inverse function was written, to get 
the normalized data (usually a large matrix) back into the original form. 

The above processes yielded several data sets that were all extremely large. These data sets were fed to an 
algorithm that built up an artificial neural net. The algorithm started with a single hidden layer and a single 
hidden neuron and attempted to train. If it didn’t train successfully within a user adjustable time frame 
(usually was set to 10000 epochs) then another neuron was added to the hidden layer and the process 
started over. There were several training methods employed also. The control was standard 
backpropagation, essentially gradient descent. The problem with this is that it is slow and prone to local 
minima in the problem space. Since the data set was nonlinear and hugely multi-dimensional and alsc > just 
huge in amount, there were a lot of chances to run into local minima. And that is just what happened. Each 
training run took several days to complete using the NASA VAB computers because of the size of the sets 
but almost never did the runs train to the minimum error criterion. Rather they stopped at the maximum 



number of hidden neurons. Each training run showed a quick drop in error initially to a plateau that just 
would not go away, no matter which training set was tried. 


The training sets were then pared down in size by about 50%. These seemed to tram well but generalization 
was a distinct problem. That is, feeding in data that the net had not seen yet (from the test set) did not yield 
the desired outputs. In fact, the output was extremely sensitive to very slight changes in input parameters 
and would output unusable controls. For example, an input set might yield a desired control of > X) degrees 
in the bank angle, a sure way to crash the simulation. 


The above was then tried with the LM training method. We were still using the same neural net architecture 
and training sets (normalized full, non-normalized full, normalized subset, non-normalized subset) but the 
training of the net was performed differently. 

LM is a method that can be considered as an improved backpropagation method. It is much faster and is 
also better at avoiding local minima in the problem space (those plateaus we ran into during^ previous 
training). But it still suffers from the heritage of being basically a gradient descent algorithm. That is a 
background on the previous work. After the above it was decided to take a step back and do an exhaustive 
examination on the computational intelligence paradigm itself. Was it the huge amount of data that was the 
problem, was it a flaw in the whole approach, or is it something else or a combination of other things . 


To this end, several strategies were examined without using the POST generated data rather more generic 
test sets provided by the literature were used. The basic architecture was examined including the activation 
functions, as were the training methods, all with a view to a more basic understanding of the process. 
Generalization, training time, and processing power were the important parameters under considera ion. 


Several small parametric studies examined parameter modification in artificial neural nets. The goals were 
to understand how these parameters affect: 


1 ) Generalization, how does the net handle inputs its never seen? 

2) Training Time, epochs needed to drive error to desired 

3) Processing Load, number of neurons needed to successfully tram 

In brief, the following areas were examined: 


1) Neural Net Architectures 

2) Test Sets 

3) Data Representation 

4) Training Methods 


Neural Net Architectures 

There are 2 basic architectures that were looked at, a single hidden layer network and a 
network. The effect of one hidden layer vs. two hidden layers was examined 
Generalization, Training Time, and Processor Load. 


dual hidden layer 
with respect to: 


1) Single hidden layer network 



2) Dual hidden layer network 



Each architecture has a varying number of nodes in each layer, above is just for Ulustrat 

The hare mtnimum of nodes for each permutation was ^ w* 

a single neuron in each hidden layer and attempting tc • train- « ^ “ h ™ dden layer would be 

added incrementally and training starte again, or e hjdd j’ er wou ld not be much larger than 

incremented first and then the first layer. This w^ J to 5000 epochs. For the single 

S^amS^ 

This 'growth' of neurons was implemented specifically “ “ J 

— — rr the - 0 ,™. numlmr of neurons ,0 

train is a big step in this. 

linear using the Matlab ‘purelin function. 

The classic activation function used for ihe hidden layer neurons is called the 'logistic sigmoid' function 
and is given by: 

/(*) = 


This function has an output between 0 and I and has a characteristic curve giving it the common name: “S- 
function”. The Matlab name for this function is logsig . 


The second function is a variation on 


the sigmoid function called the ‘hyperbolic tangent’ and looks like: 


f(x) = tanh(*) 


The function has an output between [-1.1] with the Matlab function name of 'tansig' 

The standard linear function was used for the output layers. 

f(x) = x 

Notice the difference between the logsrg and the tansig functions are only in the output range. Tansig is 
symmetric about zero while logsig is symmetric about one. 


The input/output graphs of the activation functions used look like the following figure. 



Both the sigmoid 
linear problems, 
network. 


and hyberbolic tangent are non-linear equations and are known to work well with non- 
The linear function is traditionally used in the output layer of a pattern classification 


Test Sets 

Before re-trying the POST data on the neural net, three standard test sets were examined 


1) XOR 

2) Iris 

3) Schaffer f6 


The test sets were chosen because they represent a complexity increasing standard non-linear test bed for 
computational intelligence algorithms. That is, if the methods work on these 3 known data sets then they 
should also work on the POST generated data sets or POST generated real-time data. If they don t work on 
the POST data at that point then logic would suggest that the method of gathering data should be fine tuned 
rather than the training algorithm. 


XOR 

The XOR test set is 


Input 

Output 

(0,0) ! 

0 

(0,1) 

1 

L (i.o) 

1 

d,i) 

0 


The ‘Iris’ data set is a famous set of parameters that describes 3 separate species of Iris floweis. They are 
classified by sepal width, sepal length, petal width, and petal length. There are 150 states in the set, each o 
the 3 flowers gets 50 states. So, this set can be thought of as a 4 input, single output system. The four inputs 
would be the petal/sepal parameters and the output is a number that corresponds to the flower classification 
In this case, flower 1 = 1, Bower 2 = 2, and flower 3 = 3. The Iris data set is a non-l.near mapping, hence it 
is an excellent test base for artificial neural network and more general computational intelligence research. 
The full non-nortnalized Iris Data Set is presented in the appendix. 


For example, the first 6 entries of the set in its raw form look like: 




Species 

Sepal 

length 

Sepal 

width 

Petal 

length 

Petal 

width 

1 

49 

30 

14 

2 

1 

51 

38 

19 

4 

1 

52 

41 

15 

1 

1 

54 

34 

15 

4 

1 

50 

36 

14 

2 

1 

57 

44 

15 

4 


The data was split up into three parts, representing each flower classification. Then a percentage of each of 
these three subsets was taken for training the neural nets and the rest was taken for testing. Splitting the 
full set up into 3 subsets based on each species guaranteed that each species would get equally represented 
training. 50% of the training set was used leaving 50% for testing respectively. 

Schaffer f6 

The f6 test set is derived from the f6 equation. It is highly nonlinear with a single known global minimum 
at (0,0). There are plenty of opportunities for computational intelligence algorithms to get stuck in local 
minima, so this test set is very good for testing that aspect. 

z = 0.5+ (sin 2 (sqrt(x 2 +y 2 ))-0.5)/(( i+0.01 *(x 2 +y 2 )) 2 ) 



Data Representation 

The Iris data set represents characteristics for 3 different species of flowering Iris plant. This data can be 
represented two separate ways easily in a neural network. 

First, the net can be trained to output a single number. This number will represent one of the three flowers. 
In this project the specific numbers are: flower 1 = 0.33, flower 2 = 0.66, flower 3 = 1.0. If the net output a 
number within +/-20% of the target, that was deemed close enough to consider as a match. For example, 
anything between .80 and 1.2 would equal species 3. 

Secondly, the net can be trained to output a single state. For this, three outputs were used instead of one. 
The target was: 

[1,0,0] = flower 1 
[0,1,0] = flower 2 


[0,0,1] = flower 3 

The output was considered a match if the magnitude of one of the columns was larger than the other two. 
For example, a result of [0.6,0.9,-0.1] would be interpreted as indicating flower 2 since 0.9 is larger than 
the other 2 columns. Then we just filtered for this condition with basic logic. 

In addition we can choose to normalize or not. Traditional methods call for normalizing all data The 
Matlab implementations used in the coding of this project do not require normalization in order to produce 
results. So a run was performed with normalization and without and then compared for the properties of 
generalization, training time, and processor load. 


The first six entries of normalized Iris data look like: 


0.3333 

0.3333 

0.3333 

0.3333 

0.3333 

0.3333 


0.6203 

0.6456 

0.6582 

0.6835 

0.6329 

0.7215 


0.6818 

0.8636 

0.9318 

0.7727 

0.8182 

1.0000 


0.2029 

0.2754 

0.2174 

0.2174 

0.2029 

0.2174 


0.0800 

0.1600 

0.0400 

0.1600 

0.0800 

0.1600 


It is the same data, just normalized. Similar concepts were applied to the XOR data while the f6 dataset was 
left alone since it has one minimum and is already normalized. We can relate this to the MRSR lander 
problem by recognizing that if we don't have to normalize we are better off computationally, especia y 
since normalization requires us to know the extremes if the raw data which we may "” t k "^ h f U 
normalization may make the net (hence the controller) train to new data quicker and/or with less 
computational power. It’s a trade off we need to know more about. Also we would like to know if the 
controller will work better if we had the output neurons hold a single value for bank angle or perhaps 
output neurons each with a bank angle, one representing a negative bank. 

For the normalized data sets the ‘logsig’ function was used in the hidden layers as we only need to worry 
about the range [0, 1] for the data. For the non-normalized data sets the tansig function was used in the 

hidden layers. 


Training Algorithms 

The training algorithms investigated have been: 

1) Gradient Descent, i.e. Standard Backpropagation 

2) LM (modified gradient descent) 

3) Particle Swarm Optimization (PSO) 


Standard Backpropagation algorithm, weights of the net are changed according to a difference in error 
between previous and present values. This is the control case. Prone to local minima. 


A much improved version of gradient descent. It is faster, more robust, and less computationally intensive 
overall but suffers from some of the same handicaps such as it can train any network as long as its weight, 
net input, and transfer functions have derivative functions. Less prone to local minima but still easily 

caught. 

Particle Swarm Optimization . r . . . .. . * 

In theory the PSO method is much faster and not prone to local minima. In practice of this study it is much 
slower than LM but much better at driving down to convergence though it still seems to have some 
problems with minima. An advantage of PSO is that neural trained with th.s method seem to be better at 
generalization than nets trained with traditional methods. The heart of the PSO algorithm ,s the velocity 
equation: 



Vel(k+l)=randl*vel(k)+rand2*(pbest-pos(k))+rand3*(gbest-pos(k)), k is the discrete time step 

Pos(k+ 1 )=Pos(k)+Vel(k+ 1 ) 

Pbest is the personal best valued position for the particular particle and Gbest is the global best position so 
far for entire population of particles. The variables randl, rand2, and rand3 are random numbers on the 
interval [0 1] Pos(k) is present position of a single particle, and has as many dimensions as the problem. 
This is the basic form of the PSO equation. The most common form is to include an inertia term and 
acceleration constants: 


Vel(k+1)= [inertia*rand 1 *vel(k)] 

+[ac 1 *rand2*(pbest-pos(k))] 

+[ac2*rand3*(gbest-pos(k))] 

The acceleration constants ac 1 and ac2 are fixed and have been determined by Dr. Eberhart and Dr. Shi to 
be best for most problems when both are set to 2. For our problem these values actually degrade 
performance and it seems the best values are 2 for acl and 1 for ac2. This weighs the influence ot local best 
particles more than global. There are lots of rules of thumb but not much in the way of heuristics for this 
method as it is still relatively new. That isn’t necessarily a big handicap since the algorithm is stochastic in 
nature anyway. The recommended inertia term is a linear function with time that decreases from 0.9 to 0. 
or a static term set at 1.4. Performance is definitely better with the linear term rather than the constant 
Essentially as training progresses the influence of past velocity becomes smaller. I have tried many sets ot 
the above parameters with varying degrees of success. 

For neural net training there are another couple of important parameters. The initial weight size 
(corresponds to initial positions of particles) and the maximum velocity components allowed are very 
important. Literature mentions that a maximum velocity component of 4 seems to work well but in t is 
study that is 3 orders of magnitude too large. Also, you would logically think that a large initial weight 
spread would make it easier to find the global best but experience in this research has shown the opposite is 
the case It seems that for neural net training, if we start out with a particle population tightly clustered 
together they like to travel outward to find the solution. Why this happens is an active area of investigation. 
Depending on the error over time, the maximum velocity component is modified. Sometimes during w at 
seems like a hang in a local minimum well the velocity components will saturate to the max values. It this 
happens over a series of training passes (various lengths of time have been tried) the max velocity will get a 
small decrease or increase (decrease seems to work better). This seems to jitter the problem out of local 

minima. 


Results . 

Here is the output on the Matlab workspace for a particular run (training set 50%, hidden layers as shown. 
Iris data set, normalized, training method LM): 


for 50% of the training set: 

1) Single Output/Single Hidden layer: hidden=logsig, output=purelin 
Architecture from left to right: 2 1 

Epochs to train: 33 

Generalization result: 62/75 passed test 
or 82.6667% 

2) Triple Output/Single Hidden layer: hidden=logsig, output=purelin 
Architecture from left to right: 2 3 

Epochs to train: 27 

Generalization result: 70/75 passed test 
or 93.3333% 

3) Single Output/Double Hidden layer: hidden l=purelin hidden2=logsig, output=purelin 
Architecture from left to right: 3 3 1 

Epochs to train: 49 


Generalization result: 71/75 passed test 
or 94.6667% 

4) Triple Output/Double Hidden layer: hiddenl=purelin hidden2=logsig, output=purelin 
Architecture from left to right: 6 5 3 
Epochs to train: 56 

Generalization result: 72/75 passed test 
or 96% 


All training sets followed a similar trend with virtually no difference for these data sets between using 
logsig and tansig in the hidden layers. LM outperformed Gradient Descent in all areas while PSO 
outperformed LM slightly except in time to train where it was very slow. This goes against the literature 

and warrants further study. 


V/Uliviuoiuii . J r I 

A network can have quite good results by recasting the data into multiple outputs instead of one, using only 
one hidden layer, and using either the logsig/linear or tansig/linear architectures. Overall these have the best 
performance in terms of an inductive average of the three metrics: generalization, training time, and 

processing load. 

For the MRSR problem, generalization is the most important criteria. This paper shows that a high degree 
of generalization is possible with a simple non-linear artificial neural network that uses very little 
processing power and takes almost no time to train. For any application that has the potential to injure 
people, generalization must be the overriding factor and percentages even higher than the 96% shown 
possible with this study would need to be achieved. 

If processing power and time to train (for example in a real time application) are important. Gradient 
Descent may not be a smart choice. Both LM and PSO outperformed it in those areas. 

This study confirmed that a Cl approach may be a valid and preferred approach which leads to the 
conclusion that the initial problem formulation of the MRSR along with the data preparation and mining 
need to be re-examined. 
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Appendix 

Steerable Parachute POST input deck: 
l$search 

c ****************************************************************************** 
c m2001newbase2-liftchute-targref.inp 

c 
c 

c NEW altitude target: 2700m (was 2500m), this allows relaxation of altitude 

c targeting tolerances (4-/- 150m instead of 1) and simulates about +200m above surface where we should 
c probably have a seperate terminal descent simulation anyway, 
c 
c 

c events 27 has been taken out and event 26 is a criterion 
c to stop bankangle maneuvers 
c 

c Mars Lander from LaRC hand-off to 2500 m 

c no optimization, fixed bank on steerable chute (0 deg) 

c 

c Parachute diameter = 13.0 m 
c 

c Marsgram atmosphere - Feb. 3, 2002; 00 hr, 00 min, 0.0 s 
c ln(pres) and atem input as tables 
c Marsgram winds*Braun multipliers input as tables 
c 

c This deck optimizes (min dprngl ) 
c 

c Parachute is lifting with 4 steering events 
c 1st guesses are input from bankanglestudy.dat 1 
c and are just constant turns allowing a study of achievable 
c range using this steerable chute model 
c 

c Only aoa and thrust level are controls during rev grav turn 
c 

c This deck is a modified version of EMQ-grvtn-newref.inp 
c Includes most recent M2001 chute handoff conditions & ref (south site) 
c Includes more up to date engine 

c (from CPIA/M5 Liquid propellant engine manual - unit 187): 
c Viking lander(max throttle cond): 
c Isp=210sec, Thrust=6321bf, exit area=1.588in A 2 

c Flow rate=3.10 lb/s., thrust coeff 1.53 

c 

c Constraints: 

c At event 80 (critr = wr = 2.0) - 
c 2550 < gdalt < 2850 
c -LI < ur < LI 

c -LI < vr < LI 

c 

c ******************************************************************************* 
ioflag =3, / metric input, metric output 

ipro =-l, / trajectory print flag 

c *********** *********************************************************** ********* 
c Optimization variable - 5.b 
c - general optimization inputs 

irscl =3, / default value (flag) 

isens =0, / forward finite differences, default value (flag) 

maxitr = 50, / maximum number of iteration 



opt = -i, 

optph = 80, 

optvar = ’dprngl’, 

srchm = 4, / projected gradient 

wopt =1.0, / weighting for optimization variable (default val) 

c - projected gradient specific inputs, all defaults except coneps(2) 
coneps = 89.9, 4*1.0e-04, 
consex = 1.0e-6, 0.001, 
fiterr = 1.0e-6, 0.001, 
gamax =10, 
ideb = 0, 

npad = 9.0, 4.0, 14.4494, 
pctcc = 0.3, 
pdlmax = 2.0, 
pgeps = 1.0, 
prntpd = 0, 
p2min = 1.0, 
stminp =0.1, 0.1, 
stpmax = 1.0e+10, 

^^Hej^jf:******************************************** ********************* *********** 

c Constraint variable - 5.c 
c - general Dependent Variable inputs 
depph = 80,80, 80, 80, 80, /80, 80, 
depvr = ’gdalt’/ur’/vr’, urVvrV’longVgdlaf, 
ifdeg =0,0, 0,0, 0/1,1, 
indxd = 1,2, 3,4, 5/6, 7, 
ndepv =5/7, 

c - projected gradient specific inputs 

deptl = 150.0,0.1,0.1,0.1,0.1/0.01688179,0.01688179, /target within 1km 
depval = 2700.0, 1.0, 1.0, -1.0/-1.0, 93. 8023,-15. 8384, 
idepvr = 0,1, 1,-1, -1/0,0, /constraint types 

c ****************************************^*******************=t:*****H<************ 

c Controls - 5.d 

c - general Independent Control inputs 
indph(l) = 25, 
indph(2) = 26,26, 
indph(4) = 27,27, 
indph(6) = 50,50, 
indph(8) = 28,28, 
indph( 10) = 29,29, 
c 

indvr(l) = fankpcT, 
indvr(2) = bnkpc 1 Vcritr’, 
indvr(4) = Irnkpcl’, critr’, 
indvr(6) = ’etapc 1 ’,’alppc l\ 
indvr(8) = ’bnkpcT, ’critr’, 
indvr(10) = bnkpcr, ’critr’, 
c 

nindv =11, 
pert(l) = 1.0e-4, 
pert(2) = 1.0e-4,1.0e-4, 
pert(4) = 1.0e-4,1.0e-4, 
pert(6) = 1.0e-4,1.0e-4, 
pert(8) = 1.0e-4,1.0e-4, 
pert(10) = 1.0e-4,1.0e-4, 


c 


♦include ’../../bankangleguess.daf, / just covers u(l)-u(5), i.e. chute events 
u(6) = 1. 0,0.0, 

u(8) = 5.0,5. 0, 

u(10) = -5. 0,5.0, 

c 

c - projected gradient specific input 
modew = 1, 

$ 

c Trajectory Simulation Inputs 

c ******************************************************************* ********** ** 

l$gendat 

title=0h*EMQ4iftchute-grvtn.inp*, 
event = 1 ,0.0, / first event number (primary event) 

fesn = 500, / final event number 

c - NUMERICAL INTEGRATION METHODS p. 6.a.l5-l 
npc(2) =1, / integration method (using RK) [flag] 
dltmax =1, / max step size when using variable steps [s] 
dltmin =0.05, / min step size " " " " [s] 

dt= I ,/ integration time step [s] 

kstpmx=5,/ max # of integration steps for each integration 
npinc=5,/ # of integration steps on each cycle 

^******************************************************************************* 

c Initial Event Conditions/Setup 
c 

c - INITIAL POSITION AND VELOCITY p. 6.a.l2-l 
npc(3)=l,/ initial position in xi, yi, zi 
npc(4)=U initial velocity in vxi, vyi, vzi 
npc(40)=37 reference plane for azimuth and FPA [flag] 3=? 
c 

c - initial conditions for batch runs, includes ’timeoVvia Matlab) 

^include ’../.. /pvstates.daf, 

c 

c 

c - RANGE CALCULATIONS p. 6.a. 19-1 

npc( 12)= \J cross/down range option [flag] (3=?) 

/ lonref=93. 64843 17,/ from Scott (new M2001 south site) 

/latrefgd=- 15.8 108 183, 

lonref = 93.8023, / developed from looking at averages of various 
latrefgd = -15.8384, / lifting chute non-optimized cases (showellipsebanks2.m) 
c 

c - PARACHUTE MODEL p. 6.a.28-l 
npc(32)=2, 

diamp(l)=0.0,/ init val of chute diam, unfurl to 13m 
drgpk(l)=l, 
idrgp( 1)=0, 
parif(l)=70.0, 
c 

c - AERODYNAMIC INPUTS p. 6.a.l-l 
npc(8)=3, /aerodynamic coefficient [flag] 
sref=4.5238934,/aerodynamic reference area [m A 2] (from M98) 
c 

c - AEROHEATING CALCULATIONS p. 6.a.2-l 

npc(15)=l,/ calculate aeroheating rate & tot. heat using Chapman 

npc(26)=0,/ no special aeroheating calculations 

rn= 0.6638,/ nose radius for Chapman heating (M98nom.inp) 


c - ATMOSPHERE PARAMETERS p. 6.a.4-l 

npc(5)= 6 ,/ 2/3/02, 0 hr Marsgram atmosphere input as tables 
npc(6)= 2 ,/ Marsgram winds* Braun multipliers input as tables 
atmosk(l)=241.0, 
atmosk(2)=5.335e-03, 
c 

c - CONIC CALCULATION OPTION p. 6.a.5-l 
npc(l)= 3,/Keplerian conic option [flag] 
mre=lhu, /value of mean radius to be used [m](lhu = [re+rp]/2) 
c 

c - GRAVITY MODEL p. 6.a.l0-l 

npc(16)=0,/ spherical or oblate model (oblate) [flag] 

j 2=0. 19586 16e-02,/ spherical harmonics of gravity potential function 

j 3=0.3 144926e-04, 

j4=-0. 1 889437e-04, 

j5=0.2669248e-05 , 

j6=-0.1340757e-05, 

j7=0.0d0, 

j8=0.0d0, 

re=3393940.0 ,/ equatorial radius [m] 
rp=3376780.0,/ polar radius [m] 

mu=4.28282868534e+13,/ gravitational constant (mars) [m A 3/s A 2] 
omega=7.088218e-05,/ rate of rotation of planet [rad/s] 
go=3.718,/ weight to mass factor (Mars surface) 
c - VELOCITY LOSSES p. 6.a.25-l 
npc(25)=2,/ velocity loss calculation 

c *********^***********************************************************J|C ********* 

c Initial Guidance Inputs 

^**************************************************************************** : ^** 
iguid(l)=0,0,l,/ atm.rel. aero angle guidance 
alppc(l)=180.0,/ initial alpha 
maxtim=2000.,/ maximum time 
altmax=550000.,/ maximum altitude 
altmin=-3000.0y minimum altitude 

c ********* *********************************************************** *********** 
c Vehicle Model 

C 5j£** ************************************************************* *************** 

wgtsg=2176.811,/veh. wt. at parachute deploy, N (585.479 kg) 
wpropi=372.0yinitial propellant weight, N (100 kg) 
npc(30)=07enhanced (component) weight model 
npc(9)=l, /rocket engine 

npc(27) = 1, / integrate flow rate of specified engines 
npc(22)=2, /input all four coefs in throttling parameter 
neng=2,/2 engines 

ispv(l)=553. 9, 553.97 Mars Isp (Earth Isp = 210 sec) 
iwdf(l)=2,2,/flow rate = vac. thrust/ispv 
iwpf(l)= 0,0, 

iengmf(l)=0,07engine off initially 
iengt(l)=0,07fixed engine angles (in tables) w.r.t body 
c ******************************************************************************* 
c Print Block 

c - PRINT VARIABLE REQUESTS p. 6.a.l6-l -» 
npvl=07 # of print variables per line [flag] 
pinc=10,/ print interval 
prnc=07 make profile for plotting 


include ’../../prntblk.dat’, / printing variables 
$ 

l$tblmlt 
vwum = 1.0, 
vwvm = 1.0, 

$ 

l$tab 

table = ’denkt’,1, gdalt’,3,1,1,1, 

0.0,1.0,30000,1.0,130000,1.0, 

$ 

l$tab 

table = ’prest’,1, gdalt’,6, 1,1,1, 

0,6.52553,1859,6.36383,4217,6.15610,6387,5.96265,8515,5.77084, 

8816,5.74359, 

$ 

l$tab 

table = ’atemt\l, ’gdalt’,6, 1,1,1, 

0,226. 14 1 9, 1 859,223.766,42 17,220.502,6387,2 17.854,85 15,2 15.292, 
8816,214.936, 

$ 

l$tab 

table = ’vwut’,1, ’gdalt’,6, 1,1,1, 

0,-0.1248,1859,0.02672,4217,0.19915,6387,0.55803,8515,0.86006, 

8816,0.85601, 

$ 

l$tab 

table = ’vwvt’,1, ’gdalt’,6, 1,1,1, 

0,0.200574, 1 859,0.42458,42 17,0.689 1 1 ,6387, 1 .5575,85 15,2.2659, 
8816,2.2508, 

$ 

l$tab 

table = ’vwwt’, 0,0.0, 

$ 

l$tab 

table = ’tvc It ’,0,2646.0, 

$ 

l$tab 

table = tvc2t’, 0,2646.0, 

$ 

l$tab 

table = ’ae It’, 0,0.001, 

$ 

l$tab 

table = ’ae2t’,0,0.001, 

$ 

l$tab 

table = ’pi lt’,0, 180.0, / engine #1 gimbal pitch angle 
$ 

l$tab 

table = ’yi It’, 0,0.0, / engine #1 gimbal yaw angle 

$ 

l$tab 

table = ’pi2t’, 0,0.0, / engine #2 gimbal pitch angle 

$ 

l$tab 

table = ’yi2t’,0,0.0, / engine #2 gimbal yaw angle 


$ 

l$tab 

table = ’wdlt’,0,12.98, / engine #1 flow rate 

$ 

l$tab 

table = \vd2t’, 0,12.98, / engine #2 flow rate 

$ 

l$tab 

table=’cdt’,0,1.7, 

$ 

l$tab 

table=’clt’, 0,0.0, 

$ 

l$tab 

table=’cdp 1 1’,0,0.4 1 , 
endphs=l, 

$ 

c 

c Parachute fully deployed 
l$gendat 

event=22.,0.0, / primary event 

critr= ’diamp 1 

value=13.0, 

parif(l)=0.0, 

wgtsg - 1937.297, 

endphs= 1 , 

$ 

c — Convert parachute to ’airplane’, simulates lifting/steerable chute 
l$gendat 

event=25.,0.0, /primary event 

critr=’tdurp’, 

value=0.0, 

sref= 132.73, / surface area of parachute about 13m (pi*r A 2) incorporated into lander 

npc(32)=0, / don’t calculate parachute drag 

iguid(l)=0, / aero-guidance 

iguid(2)=l, / individual component steering 

iguid(6)=0, / alpha carried over 

iguid(7)=0, / beta carried over 

iguid(8)=l, / bnkang input (from targeting algorithm) 

$ 

l$tblmlt 

$ 

l$tab 

table=’cdt ’,0,0.41, 

$ 

l$tab 

table=’clt’,0,.41, /lift coeff includes parachute, i.e. =1*0.41 
endphs = 1 , 

$ 

c - - 2nd chance to steer ’airplane’ 
l$gendat 
event=26., 
critr=’tdurp’, 
endphs = 1 , 

$ 


c - - 3rd chance to steer ’airplane’ 



l$gendat 
event=27., 
critr=’tdurp\ 
endphs = 1 , 

$ 

c - - 4th chance to steer ’airplane' 
l$gendat 
event=28., 
critr=’tdurp’, 
endphs = 1, 

$ 

c - - 5th chance to steer ’airplane’ 
l$gendat 
event=29., 
critr=’tdurp’, 
endphs = 1 , 

$ 

c jetison parachute; turn on engine #1 (start of reverse gravity turn) 

l$gendat 
event = 50,0.0, 

critr = ’gdalt’, / roving event (in case event 26 never happens) 
value = 3500.0, 

npc(32)=0, / jettison parachute 

diamp(l)=0, 

wjett = 276.702, / weight of parachute lost (N) 

sref=2.0, / new surface area 

iengmf(l) = 1,0, 
iwpf(l) = 1,0, 

iguid(l) = 0, / aero-guidance 

iguid(2) =1, / individual component steering 

iguid(6) =1, / alpha input by targeting algo 

iguid(7) = 0, / beta carried over 

iguid(8) = 0, / bnkang carried over 

$ 

l$tblmlt 

$ 

l$tab 

table=’cdt’, 0,2.0, 

$ 

l$tab 

table=’clt’,0,0.0, 

$ 

l$tab 

table=’cdp It’, 0,0.0, 
endphs = 1 , 

$ 

cc - marks 2500 meter mark (Martian surface) 
l$gendat 

event=80,0.0, / primary event - this must happen or targeting failed 

critr=’wr’, 

value = 2.0, 

iengmf(l) = 0,0, 

iwpf(l) = 0,0, 

endphs=l, 

$ 

cc 



c 

c ****** ****** *************************************************************** ***** 
c This event marks arrival at the Martian surface 
c *********** ************************************************************* ******** 

l$gendat 

event=500,0.0, / primary event 

critr=1tdurp\ 

value=0, 

endphs= 1 , 

endjob=l, 

endprb=l, 

$ 


Complete list of all files modified/used during grant cycle 
TreePrint listing of: C:\school\NASA\Docs\NASAfiles2001to2002 

C:\school\NASA\Docs\NASAfiles200Ito2002 
|| treeprint.txt 

ii 

+--dat 

| | bkb3guess.dat 
| | bkb3guess2.dat 
| | chutestates.dat 
| | chutestates2.dat 
| j fulliris.txt 
| | hoverinitcond.dat 
j j prntblk.bak 
| j prntblk.dat 
| j pvstates.dat 
1 1 

+— docs 

| | List of all papers for PHD research.doc 
j j PHD TIMELINE.doc 
| | Progress Update 1 .doc 
j j Proposalsum2001.doc 
j | Research Update.doc 
| | summer.txt 
j | Thesistemplatel.doc 

1 1 

+-Mat 

| || bc4batch.mat 
| || chutestatplot.mat 
| || copy_of_bc 1871. mat 
j II copy_of_lc.mat 
| || copy_of_lcl747_18.mat 
j || copy_of_m2001newbase.mat 
| || copy_of_rwplrps.mat 
| || forward. mat 
j jj forward3.mat 
| || inverse. mat 
j jj IrisWtsLM.mat 
j |j Ic4batch.mat 
| || lcsaveas.mat 

j jj nettrain-76init-bcl7baseline.mat 

nettrain-76init-bc!7baselineMOD.mat 



| || nettrain-nominit-bcl7baseline.mat 

| || nettrain-nominit-bcl7base!ineMOD.mat 

| || nettrainl.mat 

| || nettrain2.mat 

| || nettrain3.mat 

| || nettrainer-weights.mat 

| || nettrainer-wts-1 layer. mat 

| || nettrainer2-wts-llayer-inertialBPX.mat 

| || nettrainer2-wts-llayer-inertiallrn.mat 

| || nettrainer2-wts-llayer-inertialPS0.mat 

| || nettrainer2-wts-21ayers-unmodified-inertial-lm.mat 

| || pathdatl.mat 

| || pickstates4.mat 

| || pickstates5.mat 

| || weightsfb.mat 

| || weightsxor.mat 

i ii 

| V-oldMAT 

| | allptsCHUTEno25.mat 

| | lc4batch.mat 

i i 

+— outs 

| | copy_of_lc.out 
| | copy_of_m2001newbase.out 
| | lc4batch.out 
| | lcsaveas.out 
1 1 

| \— oldout 

| | EMQ-grvtn-newref.out 

| | hover30deg.out 

| | m2001newbase-liftchute-noopt-bnk0.out 

| | m2001newbase-liftchute-noopt-bnk25.out 

| | m2001newbase-liftchute-noopt-bnkneg25.out 

| | m2001newbase-Iiftchute-noopt-bnkstudy.out 

| | m2001newbase-Iiftchutel5LD-noopt-bnkstudy.out 

| | m2001newbase-liftchutelLD-noopt-bnkstudy.out 

| | m2001newbase-liftchute2.5LD-noopt-bnkstudy.out 

| | m2001newbase-liftchute25LD-noopt-bnkstudy.out 

| | m2001newbase-liftchute2LD-noopt-bnkstudy.out 

| | m2001newbase-liftchute75LD-noopt-bnkstudy.out 

| | m2001newbase-thrustchute-noopt-bnkstudy.out 

| | m200inewbase-thrustchutel2kg-noopt-bnkstudy.out 

| | m2001newbase.out 

| | m2001newbase2-liftchute-noopt-bnkstudy.out 

| | m2001newbase2-liftchute-targref.out 

| | m2001newbase2-liftchute.out 

| | m2001newbase2-liftchutel5LD-noopt-bnkstudy.out 

| | m2001newbase2-liftchutel5LD-targref.out 

j j m2001newbase2-liftchutel5LD.out 

| | m2001newbase2-liftchutelLD-noopt-bnkstudy.out 

| | m2001newbase2-liftchutelLD-targref.out 

| j m2001newbase2-liftchutelLD.out 

| | m2001newbase2-liftchute25LD-noopt-bnkstudy.out 

| j m2001newbase2-liftchute25LD-targref.out 

j | m200 1 ne wbase2-liftchute25LD.out 

j | m2001newbase2-liftchute2LD-noopt-bnkstudy.out 



| | m2001newbase2-liftchute2LD-targref.out 

j j m2001newbase2-liftchute2LD.out 

| | m200inewbase2-liftchute75LD-noopt-bnkstudy.out 

j j m2001newbase2-liftchute75LD-targref.out 

j j m2001newbase2-liftchute75LD.out 

| | m200 1 ne wbase2-thrustchute-noopt-bnkstudy .out 

| | m2001newbase2-thrustchute.out 

] | m2001newbase2-thrustchutel2kg~noopt-bnkstudy.out 

| | m2001newbase2-thrustchute 12kg. out 

| | m2001newbase2.out 

| | tstnet.out 

i i 

+— post3d_source 
| aero.f 
| aero.o 
| balland.f 
| balland.o 
| calspe.f 
| calspe.o 
| guidance, f 

| guidance, o 

| Makefile 
| Makefile.log 
| Makefile_calspe 
| Makefile_calspe_new 
| Makefile_calspe_oId 
j Makefile.dbg 
| Makefile_dbg„calspe 
| mars_entry5a.inp 
| master, o 
I P3d 

| p3d_calspe 
| p3d_dbg 
| p3d_dbg_calspe 
| phzxm.o 

i 

V-support 

chutestatplot.m 

demonormalize.m 

extho.m 

f6.m 

f6pso.m 

learnga.m 

learnlm.m 

load_D AT A_mat. m 

m2001viewsite.m 

MatPathData.m 

nettrainer.m 

nettrainer2.m 

normalize, m 

p2check.m 

p3dbatchSPACE2000.m 

p3dbatchspace2000PSO.m 

p3dSPACE2000.m 

p3dSPACE2000b.m 

pso.m 


startup. m 

tbpxl.m 

tbpx2.m 

tbpx3.m 

testf6.m 

testnettrainer.m 

testpso.m 

tgal.m 

tga2.m 

tlml.m 

tlm2.m 

tlm3.m 

tpsol.m 

tpso2.m 

tpso2mod.m 

tpso2old.m 

tpso3.m 

trainbpx.m 

trainga.m 

trainlm.m 

trainpso.m 

tstmov.m 

unwrapmat.m 

wrapmat.m 


Iris Data Set 

The Iris data set used here was downloaded from a website after a google search for the words “Iris Data 
Set'. It is presented here in the exact downloaded form: 
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